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Preface 


This  report  presents  the  results  of  a  numerical  model  investigation  used  to 
predict  average  salinity  changes  that  will  occur  in  the  Galveston  Bay  region  as  a 
result  of  the  addition  of  barge  lanes  along  the  Houston  Ship  Channel. 

This  investigation  was  conducted  from  September  2000  through  May  2001  at 
the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  by  Dr.  R.  C. 
Berger,  Mr.  A.  R.  Carrillo,  Ms.  M.  S.  Sarruff,  and  Mr.  J.  A.  Sanchez  of  the 
Coastal  and  Hydraulics  Laboratory  (CHL).  Funding  was  provided  by  the  U.S. 
Army  Engineer  District,  Galveston. 

The  work  was  performed  under  the  general  direction  of  Mr.  Thomas  W. 
Richardson,  Director,  CHL,  and  Dr.  R.  T.  McAdory,  Chief,  Tidal  Hydraulics 
Branch,  CHL. 

At  the  time  of  publication  of  this  report,  Dr.  James  R.  Houston  was  Director 
of  ERDC,  and  COL  John  W.  Morris  IE,  EN,  was  Commander  and  Executive 
Director. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  of  such  commercial 
products. 


Conversion  Factors,  Non  SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  meters 

feet 

0.3048 

meters 

miles  (U.S.  statute) 

1.609347 

kilometers 

square  feet 

0.09290304 

square  meters 

1 


Introduction 


Background 

The  U.S.  Army  Engineer  District,  Galveston,  is  making  an  investigation  of 
the  addition  of  barge  lanes  along  the  Houston  Ship  Channel  through  Galveston 
Bay,  TX.  In  the  early  to  mid  1990s  the  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC),  Coastal  and  Hydraulics  Laboratory  (CHL) 
produced  a  three-dimensional  (3-D)  numerical  model  of  Galveston  Bay  to  study 
the  enlargement  of  the  Houston  Ship  Channel  from  nominal  dimensions  of  40  ft1 
deep  mlw  by  a  400-ft-width  (base)  to  45-ft-deep  channel  with  a  530-ft  width 
(phase  I)  (Berger  et  al.  1995a).  The  model  produced  water  surface,  velocity,  and 
salinity  fields  for  the  base  and  phase  I  conditions.  These  data  in  turn  were  used  to 
make  an  environmental  assessment.  The  channel  is  presently  being  dredged  to 
phase  I  dimensions. 

The  proposed  barge  lanes  are  to  allow  barge  traffic  to  safely  coexist  with  the 
deep-draft  navigation.  These  lanes  are  to  be  along  and  on  both  sides  of  the 
channel  throughout  Galveston  Bay.  The  barge  lanes  are  to  have  a  1 4-ft  required 
depth  and  extend  500  ft  from  the  center  line  of  the  navigation  channel.  The 
Galveston  District  requested  that  the  Galveston  Bay  model  be  used  to  provide 
salinity,  velocity,  and  water  surfaces  in  support  of  an  environmental  assessment. 

The  Galveston  Bay  model  consists  of  the  finite  element  mesh,  the  data,  and 
the  hydrodynamic  code.  The  finite  element  mesh  represents  the  spatial  grid  and 
bathymetry.  The  model  was  updated  to  match  the  “as  built”  conditions.  These 
include  the  actual  disposal  configuration,  the  Bayport  Channel  bathymetry,  and 
the  actual  dredged  channel  depths  for  the  Houston  Ship  Channel.  The  data 
includes  the  boundary  conditions  and  forcing  data  (winds,  Coriolis,  etc.).  The 
hydrodynamic  code  used  has  evolved  since  the  prior  investigation  with  the 
addition  of  sediment  transport,  Smagorinsky  horizontal  diffusion,  and  improved 
Newton  iteration,  among  several  items.  That  particular  code  is  now  referred  to  as 
TABS-MDS. 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  vi. 
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Objective 


The  objective  of  the  work  presented  herein  is  to  predict  the  month-average 
salinity  changes  that  will  occur  in  Galveston  Bay  (Figure  1 ),  resulting  from  the 
addition  of  the  proposed  barge  lanes.  The  tasks  undertaken,  and  presented  herein 
are: 

a.  Provide  a  short  check  of  the  verification  using  the  TABS-MDS  code. 

b.  Incorporate  all  features  of  the  as  built  phase  I. 

c.  Modify  the  grid  to  include  the  barge  lane  configuration. 

d.  Run  three  plans  for  one  hydrologic  year  each.  The  three  plans  are: 

(1)  Phase  I  as  studied  in  the  prior  model  study. 

(2)  Phase  I  as  built. 

(3)  Phase  I  as  built  with  the  addition  of  the  barge  lanes. 

Since  a  low-flow  hydrologic  scenario  is  regarded  as  the  most  critical,  the 
“present”'  low-flow  boundary  conditions  are  utilized. 

The  purpose  of  this  report  is  to  present  the  results  of  a  numerical  model 
investigation  addressing  this  objective.  Galveston  Bay  is  the  modeled  region 
including  Trinity  Bay,  East  Bay,  and  West  Bay,  as  shown  in  Figure  1 . 


Figure  1.  Vicinity  map 
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Model 


Numerical  Model 

TABS-MDS,  the  numerical  model  used  for  this  investigation,  is  a  3-D  finite 
element  code  originally  developed  by  Dr.  Ian  King  of  Resource  Management 
Associates,  and  modified  at  ERDC.  It  models  3-D  hydrodynamics  and  salt 
transport  accounting  for  unsteady  river  inflows,  tides,  wind  effects,  and  density- 
driven  circulation.  It  has  been  widely  used  by  ERDC  to  model  3-D 
hydrodynamics  and  salinity  at  numerous  locations,  including  Galveston  Bay,  TX 
(Berger  et  al.  1995a,  1995b).  The  model  is  described  more  fully  in  Appendix  A. 

Verification 

The  TABS-MDS  numerical  model  has  evolved  since  the  prior  Galveston  Bay 
verification,  so  a  short  reverification  check  was  performed  using  the  same 
boundary  conditions  and  bathymetry  as  the  original  study  (Berger  et  al.  1995a). 
This  will  provide  the  District  with  an  up-to-date  model  that  can  be  used  for  any 
hydrodynamic  modeling  of  the  bay.  A  brief  discussion  of  the  specifics  of  the 
experimental  conditions  are  addressed  in  the  following  sections.  For  a  more 
detailed  discussion  see  Berger  et  al.  1995a. 

Computational  mesh 

The  verification  check  consists  of  using  the  current  updated  TABS-MDS 
code  for  the  verification  period  used  in  the  original  model  verification  (Berger  et 
al.  1995a).  The  channel  configuration  for  the  verification  check  is  that  of  the  40 
x  400-ft  channel  that  existed  when  the  field  data  were  collected.  Figure  2 
illustrates  the  planform  view  of  the  computational  mesh  used.  The  mesh 
represents  the  1990  existing  conditions.  The  grey  areas  denote  the  3-D  regions  of 
the  mesh.  The  mesh  consists  of  12,710  nodes  and  4,422  elements.  (The  mesh 
was  refined  near  the  channel  compared  to  the  earlier  model  study.) 
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Figure  2.  Computational  mesh  -  1 990  existing  conditions  for  verification 


Verification  field  data 

Field  data  collection  occurred  from  19  July  1990  to  15  January  1991.  The 
bulk  of  these  data  was  collected  by  moored  salinity  and  velocity  meters,  as  well  as 
tide  gauges  mounted  throughout  the  system  (Figures  3  and  4).  Within  this  6- 
month  time,  a  short-duration  intense  field  data  collection  was  conducted  19-20 
July  1990.  These  data  included  velocity  and  salinity  measurements  over  depth  at 
five  ranges  along  the  navigation  channel  (See  Figure  5).  Complete  details  of  the 
field  collection  may  be  found  in  Fagerburg  et  al.  (1994). 
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Figure  3.  Basic  long-term  survey  data  collection  equipment  locations  (sta  SI  .0, 
S2.0,  S5.5,  S6.0,  S8.1,  S9.0,  S1 1.1,  SI  2.0,  and  SI  5.1  record  salinity 
as  well) 
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Figure  4.  Supplementary  long-term  data  collection  equipment  locations 
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Experimental  conditions 

The  data  used  to  drive  the  model  tide  were  the  National  Ocean  Service/ 
National  Oceanic  and  Atmospheric  Administration  (NOS/NOAA)  sta  877-1510. 
Pleasure  Pier  tide  gauge  data  shifted  forward  1.31  hr  to  account  for  the  model 
boundary  approximately  26  miles  offshore  from  Pleasure  Pier.  The  actual  data 
were  filtered  to  remove  signals  of  periods  less  than  3  hr.  The  salinity  at  the 
boundary  was  estimated  by  using  the  published  averages  supplied  by  Cochrane 
and  Kelly  (1986).  The  wind  data  were  provided  by  two  sources.  ERDC 
established  a  meteorological  station  at  location  SI 0.1  (Figure  4)  for  period  19 
July  1990  to  15  January  1991.  The  National  Weather  Service  data  at  Houston 
International  Airport  were  also  obtained  for  this  period.  A  correlation  was 
derived  between  the  National  Weather  Service  and  ERDC  data  so  that  readings  at 
Houston  International  could  be  used  for  Galveston  Bay  (see  Berger  et  al.  1995a). 
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The  model  must  be  run  for  a  length  of  time,  termed  spin-up,  to  remove  initial 
conditions  error.  The  closer  that  the  initial  salinity  field  is  to  the  ultimate  field, 
the  less  length  of  time  is  required  for  spin-up.  Therefore,  the  verification  initial 
salinity  field  for  January  1990  was  generated  by  interpolating  actual  field  data. 
This  rough  approximation  was  sufficient  since  the  verification  field  data  occurred 
in  July,  giving  a  sufficient  relaxation  time  for  the  data. 

The  numerical  model  calculated  water-surface  elevations,  current  velocities 
(3-D  components),  and  salinities  at  each  node  every  30  min  for  the  1 2-month 
period  of  simulation  between  January  1  and  December  31. 


Harmonic  analysis 

A  tidal  signal  can  be  affected  by  a  number  of  factors,  including  astronomical 
as  well  as  meteorological  ones.  Tidal  behavior  is  therefore  better  understood 
after  removing  variable  unpredictable  events,  such  as  wind-driven  magnitudes 
and  random  meteorological  effects.  A  harmonic  analysis  of  the  tidal  record  can 
remove  the  factors  mentioned,  while  exposing  the  major  tidal  components.  The 
resulting  amplitude  and  phase  prove  to  be  useful  tools  when  evaluating  the 
performance  of  a  numerical  model,  related  to  tide  approximations. 

Sta  S1.0,  S3.0,  and  SI 6.0  were  selected  for  analysis  (See  Figure  3).  An  8- 
week  period  spanning  August  25  through  October  27,  1990,  was  extracted  from 
each  station  record.  The  model  data  were  in  0.5-hr  increments,  while  the 
prototype  data  varied.  That  is,  sta  S  1.0  was  recorded  every  0.5  hr;  sta  S3.0  was 
recorded  in  0.24-hr  increments.  Finally,  sta  S 1 6.0  was  recorded  in  0.25-hr 
increments.  The  data  were  filtered  to  remove  periods  less  than  3  hr  and  greater 
than  35  hr.  The  filtering  process  yielded  data  sets  comprising  primarily  diurnal 
and  semidiurnal  components,  as  well  as  overtides. 

The  harmonic  analysis  results  in  a  list  of  amplitudes  and  phases 
corresponding  to  diumal  and  semidiurnal  constituents.  For  the  purposes  of  this 
study,  only  diumal  constituents  were  used  in  our  comparison  of  prototype  and 
model.  After  observing  the  resulting  constituents,  two  diumal  components 
dominated  the  tidal  range,  K1  and  01,  corresponding  to  periods  of  23.93  and 
25.92  hr,  respectively.  A  ratio  of  the  amplitude  of  K1  and  01  at  each  station  to 
that  of  sta  SI  .0,  or  relative  amplitude,  was  calculated  and  charted  on  Figures  6 
and  7. 
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Figure  6.  Prototype/model  tidal  amplitude  comparison  for  K1  constituent 


Figure  7.  Prototype/model  tidal  amplitude  comparison  for  Ol  constituent 


For  a  period  of  23.93  hr  (Kl),  a  significant  reduction  of  the  tidal  range  was 
evidenced  from  sta  S1.0  to  S3.0  (8.1  miles),  about  27  percent.  Upstream  to 
sta  SI 6.0  (26.4  miles),  the  tidal  amplitude  remained  steady,  increasing  slightly  by 
2  percent.  The  model  followed  this  trend  remarkably  well,  separated  from  the 
prototype  amplitude  by  only  12  percent  at  the  most.  Likewise,  the  reduction  of 
the  tidal  range  from  sta  S 1 .0  to  S3.0  was  about  35  percent  for  a  25.82  hr  period 
(01),  slightly  larger  than  that  of  constituent  Kl .  The  tidal  range  then  increased 
by  9  percent  to  sta  SI 6.0.  In  this  case,  the  model  was  even  closer  to  the  prototype 
amplitude,  separated  by  only  5  percent  at  the  most. 
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Additionally,  comparing  the  phase  lag  provided  further  insight  about  the 
model  approximations.  The  phase  lag  is  the  difference  between  the  phase  at  each 
station  and  that  of  sta  SI  .0.  Figures  8  and  9  depict  this  information  graphically. 


Figure  8.  Comparison  of  prototype  and  model  phase  lag  for  K1  constituent 


Figure  9.  Comparison  of  prototype  and  model  phase  lag  for  01  constituent 


At  sta  S3.0,  the  prototype  led  the  model  by  47  min  for  constituent  K1 . 
Upstream  at  sta  SI 6.0,  the  prototype  led  the  model  by  much  less,  only  around 
7  min.  For  constituent  01 ,  the  prototype  led  the  model  by  about  33  min,  while  at 
26.4  miles  upstream  (SI 6.0),  the  model  led  the  prototype  by  6  min. 
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Salinity 


Plate  1  shows  the  modeled  salinities  versus  the  field  data  for  three 
representative  meter  locations  (sta  6.1,  sta  12.1,  and  sta  14.0).  These  three 
stations  are  shown  as  a  short  check  example  to  demonstrate  that  the  model  is 
functioning  much  as  it  did  in  the  previous  study.  The  period  of  time  represented 
is  19  July  1990  to  15  January  1991.  A  high  flow  event  occurred  in  the  Trinity 
Basin  in  June,  so  the  data  include  significant  salinity  variation  as  the  system 
rebounded  from  near  fresh  conditions.  The  continuous  meter  readings  from  the 
field  are  shown  as  red  lines.  These  meters  were  serviced  periodically,  at  which 
time  a  single  grab  sample  was  taken.  These  are  shown  as  an  X  for  each  recorded 
value.  The  model  tracks  the  rebound  in  each  case.  For  stas  12.1  and  14.0,  the 
model  underpredicts  the  salinity  by  about  3  ppt,  whereas  at  sta  6. 1  the  rebound 
value  of  salinity  is  quite  close. 


Conclusions 

The  tidal  and  salinity  results  for  the  verification  compared  well  to  the 
previous  study.  The  model  thus  remains  calibrated  and  suitable  for  salinity  and 
hydrodynamic  studies  of  Galveston  Bay,  particularly  base-plan  type  comparisons 
as  performed  in  the  current  study. 
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3  Approach 


Plans 

Since  the  completion  of  the  channel  deepening  study,  the  actual  deepening 
operation  has  been  undertaken.  The  actual  project  has  some  differences  from 
what  was  originally  tested.  This  included  the  arrangement  of  the  disposal  islands 
as  well  as  the  actual  depth  dredged  in  the  channel.  For  this  reason,  three 
configurations  of  the  Houston  Ship  Channel  through  Galveston  Bay  were 
modeled.  The  tests  were  as  follows: 

a.  “Final”  -  phase  I  as  studied  in  the  prior  model  study. 

b.  “Built”  -  phase  I  as  built. 

c.  “Bargelanes”  -  phase  I  as  built  with  the  addition  of  the  barge  lanes. 

The  “Final”  plan  uses  the  phase  I  conditions  as  studied  in  the  prior  model 
investigation  (Berger  et  al.  1 995b),  and  has  a  47-ft-  (45  ft  authorized  plus  2  ft  of 
allowable  overdepth)  deep  channel  with  a  530-ft  width.  The  “Built”  plan  has  a 
47-  to  49-ft  deep  channel  as  was  actually  dredged.  The  “Bargelanes”  plan  adds 
the  barge  lanes  (14-ft  depths  extending  500  ft  from  the  center  line  of  the 
navigation  channel)  to  the  “Built”  plan. 

The  numerical  model  calculates  water-surface  elevations,  current  velocities 
(3-D  components),  and  salinities  at  each  node  every  30  min  for  the  1 2-month 
period  of  simulation  between  January  1  and  January  3 1 .  Those  data  are 
processed  to  provide  average  monthly  salinity  contour  plots  for  each  plan  for  each 
of  the  12  months  (Plates  2-25).  For  comparison  purposes,  the  “Built”  and 
“Bargelanes”  plans  are  plotted  on  the  same  graph. 

Computational  Meshes 

Figures  10  through  12  illustrate  the  planform  view  of  the  computational 
meshes  used.  The  grey  areas  denote  the  3-D  regions  of  the  meshes.  Table  1 
shows  the  number  of  nodes  and  elements  for  each  mesh. 
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San  Jacinto  River 


Figure  1 1 .  Computational  mesh  -  phase  I  as  built 
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San  Jacinto  River 


Figure  12.  Computational  mesh  -  phase  I  as  built  with  addition  of  barge  lanes 
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Experimental  Conditions 


Since  most  of  this  work  is  based  on  prior  investigations,  only  a  brief 
discussion  of  the  specifics  of  the  experimental  conditions  are  addressed  in  the 
following  sections.  For  a  more  detailed  discussion  see  Berger  et  al.  1995b. 


Inflows 

Modeling  conditions  for  the  plan  experiments  use  the  “Present”  low-flow 
hydrologic  era  from  the  prior  investigation.  The  low-flow  scenario  was  chosen 
since  it  is  regarded  as  the  most  critical. 

In  addition  to  the  river  inflows,  two  Reliant  Energy  electric  power  generation 
plants  were  included  in  the  model:  the  Cedar  Bayou  and  Robinson  plants. 
Monthly  intake  rates  for  the  year  1990  were  used  for  all  testing.  Salinity  values 
obtained  at  the  intake  are  used  as  the  boundary  condition  at  the  discharge. 


Tides 

The  tidal  boundary  data  were  from  NOS/NOAA  sta  877-1510  at  Pleasure 
Pier  for  calendar  year  1984,  which  was  found  to  be  a  typical  temperature  year. 
The  model  boundary,  which  is  about  26  miles  offshore,  used  these  data  shifted 
forward  in  time  by  1.31  hr.  This  matched  the  phase  at  Pleasure  Pier.  Periods  of 
less  than  3  hr  were  removed  by  filtering,  preventing  noise  in  the  results. 


Wind 

Wind  data  for  1 984  were  obtained  from  the  National  Weather  Service  for 
Houston  Intercontinental  Airport.  From  a  correlation  with  data  collected  at 
sta  SI  0.1  in  Galveston  Bay  during  the  verification  period,  these  wind  data  were 
transposed  to  the  bay  system. 


Salinity 

Gulf  boundary  salinities  are  based  on  15-year  monthly  averages  for  the  Gulf, 
and  the  work  of  Conchrane  and  Kelly  (1986).  The  final  salinity  field  results  from 
the  verification  run  are  used  as  the  initial  conditions  for  each  of  the  plan 
experiments. 
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4  Results 


Phase  I  Final  (Plates  2-13) 

This  is  the  bathymetry  and  disposal  configuration  that  matches  that 
previously  tested  in  the  Houston-Galveston  Ship  Channel  Deepening  Study 
(Berger  et  al.  1995b).  The  plates  represent  the  bottom  salinity  concentrations. 
Beginning  in  January  the  isohalines  generally  show  increasing  intrusion  along  the 
channel  until  April.  The  freshwater  inflow  reaches  its  peak  in  April  and  May  and 
the  salinity  recedes  a  bit.  The  impact  of  the  Bayport  Ship  Channel  is  sometimes 
evident  in  the  shape  of  the  contours  at  certain  times,  such  as  February,  May,  and 
June.  June  appears  to  be  the  minimum  month-average  salinity.  The  intrusion 
increases  as  the  fresh  water  is  reduced  throughout  the  summer  and  fall  months. 
The  peak  intmsion  is  in  October  and  November. 

Phase  I  as  Built  and  Phase  I  as  Built  with  Barge 
Lanes  (Plates  14-25) 

Generally  the  intrusion  follows  the  same  pattern  as  that  found  in  “Phase  I 
Final.”  The  minimum  intrusion  is  in  June  and  the  maximum  is  in  October.  The 
barge  lanes  are  14  ft  deep  at  mlw.  The  southern  reach  of  the  bay  below  Eagle 
Point  is  deeper  than  this  and  so  no  deepening  takes  place  there.  Above  Eagle 
Point  the  barge  lanes  must  be  dredged.  The  most  noticeable  differences  between 
the  conditions  with  and  without  the  barge  lanes  occurs  in  October,  November, 
and  December.  As  with  the  channel  deepening  the  primary  area  of  influence  is 
the  upper  western  portion  of  the  bay  from  Clear  Lake  and  above.  The  isohalines 
are  also  shifted  upstream  in  Trinity  Bay  in  October  and  November.  This  is 
partially  due  to  the  increased  salinity  at  the  intake  of  the  Reliant  Energy  electric 
generation  plant  at  Cedar  Bayou,  that  has  its  discharge  in  the  upper  (northern) 
Trinity  Bay,  and  to  greater  intmsion  along  the  Houston  Ship  Channel. 
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Table  2 

Modeled  Month  Average  Salinities,  for  Phase  I  as  Built  and  Phase  I  as  Built  with  Barge 
Lanes,  for  Select  Grid  Locations,  ppt 


Location  3 


Location  4 


Jan. 

Feb. 

Mar. 

Apr. 

May 

June 

July 

Aug. 

Sept. 

Oct. 

Nov. 

Dec. 

15.8 

13.8 

15.2 

15.4 

13.6 

14.2 

15.8 

18.3 

19.7 

19.8 

19.3 

17.7 

15.8 

13.7 

14.0 

15.8 

14.3 

14.3 

16.5 

18.7 

20.1 

20.7 

20.4 

18.7 

10.9 

9.0 

10.8 

10.8 

10.6 

12.6 

14.7 

15.4 

16.3 

15.4 

12.5 

11.0 

B.8 

10.2 

11.3 

10.3 

10.7 

12.9 

14.8 

15.8 

17.1 

16.1 

13.5 

8.6 

7.3 

7.1 

7.1 

EM 

6.8 

12.4 

12.9 

14.0 

14.1 

10.3 

8.5 

8.1 

7.0 

7.8 

6.6 

7.5 

nm 

11.7 

13.4 

15.1 

14.8 

10.9 

12.2 

10.6 

10.3 

10.9 

9.3 

9.8 

11.5 

14.2 

15.8 

16.6 

16.5 

14.9 

12.1 

10.8 

10.4 

10.6 

10.2 

10.5 

11.6 

14.3 

16.0 

17.3 

17.6 

15.3 

The  isohalines  shown  in  the  plates  are  contours  of  equal  concentration. 

These  are  good  for  getting  overall  salinity  distributions  and  shifts  in  salinity 
patterns,  but  can  be  difficult  in  determining  salinity  at  locations.  Table  2  shows 
the  actual  mean  monthly  values  of  salinity  at  four  locations  for  all  months  and  for 
Phase  I  as  Built”  and  “Phase  I  as  Built  with  Barge  Lanes.”  These  four  locations 
are  identified  in  Figure  13.  Locations  1  and  2  are  about  the  same  distance  up  the 
navigation  channel,  with  location  1  west  of  the  channel  and  location  2  east. 
Location  2  is  close  to  the  intake  of  Reliant  Energy's  Cedar  Bayou  Electric 
Generating  Station.  Location  3  is  further  up  Trinity  Bay  closer  to  the  discharge 
from  the  power  station.  Location  4  is  in  the  middle  of  Trinity  Bay  and  should  be 
useful  in  determining  any  additional  intrusion  directly  from  Houston  Ship 
Channel  with  the  proposed  barge  lanes. 
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Figure  13.  Select  grid  locations  for  comparison  of  phase  I  as  built  and  phase  I  as 
built  with  barge  lanes 


First  considering  location  4,  the  maximum  increase  in  month-average  salinity 
due  to  the  barge  lanes  is  1 . 1  ppt  in  November.  Location  1  shows  decrease  in 
salinity  due  to  the  barge  lanes  in  the  early  spring,  but  by  November  it  also  shows 
a  1.1  ppt  increase  in  salinity  due  to  the  barge  lanes.  Location  2  shows  a  slight 
decrease  in  salinity  in  early  spring  (0.6  ppt)  and  a  maximum  increase  of  1.0  ppt  in 
December.  Location  3  should  reflect  the  influence  of  the  outflow  from  the  power 
plant  as  well  as  the  direct  change  in  the  salinity  from  the  navigation  channel  (as 
indicated  by  location  4).  The  maximum  increase  in  salinity  is  1.1  ppt  in  October. 
This  is  closer  to  the  timing  of  the  increase  at  location  4,  which  may  suggest  that 
the  direct  intrusion  up  Trinity  Bay  is  more  influential  than  the  power  plant. 
However,  the  plates  of  isohalines  sometimes  show  the  distinctive  flow  of  salt 
water  into  upper  Trinity  Bay  from  the  power  plant  (see  Plate  23). 

The  average  change  in  salinity  is  +0.4,  +0.3,  +0.4  and  +0.4  ppt  for  locations 
1,  2,  3,  and  4  respectively.  The  most  noteworthy  increases  are  late  in  the  year  just 
after  the  peak  salinity  period.  To  get  an  idea  on  the  relative  magnitude  of  this 
change,  we  will  compare  these  results  to  that  caused  by  the  river  inflows  to  the 
bay.  From  Berger  et  al.  1995a,  a  salinity  change  of  about  1 .0  ppt  in  upper  Trinity 
Bay  can  be  developed  from  a  10  percent  change  in  freshwater  inflow.  These  0.4 
ppt  increases  in  salinity  on  average  have  roughly  the  impact  of  a  4  percent 
reduction  in  freshwater  inflow. 
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5  Summary  and  Conclusions 


The  basic  goal  of  this  modeling  effort  is  to  evaluate  the  impact  of  the  addition 
of  barge  lanes  on  each  side  of  the  Houston  Ship  Channel  through  Galveston  Bay. 

The  parameter  of  concern  is  salinity.  The  following  steps  were  taken  to  perform 
this  evaluation: 

a.  The  previously  developed  RMA10-WES  3D  hydrodynamic  and  salinity 
model  used  to  evaluate  the  Houston-Galveston  Ship  Channel  deepening 
would  be  activated  and  upgraded  to  the  present  version  of  RMA10-WES 
called  TABS-MDS. 

b.  The  mesh  was  refined  along  the  channel  to  more  accurately  represent  the 
region  in  which  the  barge  lanes  would  be  added. 

c.  The  model  was  run  through  a  verification  check,  by  using  the  1 990 
bathymetric  conditions  and  running  the  wind,  freshwater  inflow  and  tide 
conditions  used  for  the  prior  model  verification. 

d.  All  testing  after  verification  was  conducted  for  what  are  believed  to  be 
critical  conditions.  This  is  the  hydrologic  scenario  termed  1990  low 
flow. 

e.  A  single  yearlong  run  was  made  of  what  is  called  the  “Phase  I  Final” 
condition.  This  is  the  bathymetry  and  disposal  island  configuration  that 
was  run  as  the  final  condition  in  the  prior  channel  deepening  study.  The 
channel  dimensions  were  47  ft  mlw  by  530  ft  wide.  These  results  are  to 
be  used  to  make  sure  that  the  results  are  similar  to  that  given  during  the 
previous  study. 

/  Two  runs  were  then  made  over  the  same  yearlong  period.  These  were  the 
“Phase  I  as  Built”  and  the  “Phase  I  as  Built  with  Barge  Lanes  Condition.” 
The  main  navigation  channel  dimensions  were  530  ft  wide,  but  the  depth 
was  47  to  49  ft  deep  at  mlw  in  the  lower  bay.  The  upper  bay  and  the 
bayou  section  had  a  depth  of  47  ft  mlw.  The  barge  lane  depth  was  14  ft 
mlw.  These  runs  are  the  basis  of  the  salinity  evaluation  of  the  barge 
lanes. 

The  model  tide  and  salinity  verification  check  showed  that  the  present  model 
is  performing  with  precision  similar  to  that  of  the  prior  study.  The  “Phase  I  Final” 
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test  showed  that  the  present  model  predicts  a  salinity  regimen  much  like  that  of 
the  previous  study. 

The  addition  of  the  barge  lanes  had  little  impact  on  the  lower  (southern) 
portion  of  Galveston  Bay.  The  upper  region  of  Galveston  Bay  and  Trinity  Bay  on 
average  showed  about  a  0.4  ppt  increase.  In  these  same  upper  regions  the  spring 
period,  in  which  the  salinity  rebounds,  showed  a  drop  in  salinity  compared  to  the 
same  channel  configuration  without  the  barge  lanes. 

The  addition  of  the  barge  lanes  tended  to  raise  the  salinity  in  the  upper 
western  portion  of  Galveston  Bay  and  in  Trinity  Bay  in  late  fall  (October, 
November,  and  December).  This  Trinity  Bay  increase  appears  to  be  due  to  the 
increased  intmsion  directly  from  the  navigation  channel  up  Trinity  Bay,  as  well  as 
the  transfer  of  increased  salinity  at  the  intake  of  the  Reliant  Energy’s  Cedar  Bayou 
Electric  Generation  station.  This  plant  has  an  intake  near  Cedar  Bayou,  which  is 
relatively  near  the  navigation  channel  but  has  a  discharge  in  upper  Trinity  Bay. 
The  change  in  salinity  at  the  intake  in  turn  changes  the  salinity  discharged  into 
Trinity  Bay. 
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Appendix  A 
Hydrodynamic  Code 


The  complexity  of  this  estuary  requires  a  numerical  model  that  relies  upon  an 
unstructured  computational  mesh.  The  code  chosen  is  the  Galerkin-based  finite 
element  model  TABS-MDS,  which  is  a  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC)  adaptation  of  the  RMA-10  code  developed  by  King 
(1993).1  This  code  computes  time- varying  open-channel  flow  and  salinity/ 
temperature  transport  in  one,  two,  and  three  dimensions.  It  invokes  the 
hydrostatic  pressure  and  mild  slope  assumption.  Vertical  turbulence  is  supplied 
using  a  Mellor-Yamada  Level  II  (Mellor  and  Yamada  1982)  k-l  approach  modi¬ 
fied  for  stratification  by  the  method  of  Henderson-Sellers  (1984).  The  salin¬ 
ity/density  relationship  is  based  upon  Pritchard  (1982). 

The  full  3-D  equations  are  reduced  to  a  set  of  two  momentum  equations,  an 
integrated  continuity  equation,  a  convection-diffusion  equation,  and  an  equation 
of  state.  The  simplification  is  a  result  of  the  hydrostatic  pressure  approximation. 
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1  All  references  cited  in  this  appendix  are  listed  in  the  References  section  at  the  end  of  the 
main  text. 
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P=F(s) 


(A5) 


Elevation-related  terms  are  defined  in  Figure  Al. 


Figure  Al .  Definitions  for  elevation  terms 


where 


p  =  density 

u,v,w  =  x,y,z  velocity  components 
t  =  time 
P  =  pressure 

,2.2 

f>  +  +rH,>cos(ej 

c 2 

,2.2  \<!/2> 

ry=  -f£iu-pgVa(Ua+Va)  +v,w2  sinr©; 

c2 

Q  =  2(osin(^) 

co  =  rate  of  angular  rotation  of  the  earth 
cp  =  local  latitude 
g  =  gravitational  acceleration 
C  =  Chezy  or  Manning  friction  formulation 
=  a  coefficient  from  Wu  (1980) 

W  =  wind  speed 

0  =  wind  direction  counterclockwise  from  easterly 
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h  = 
uc,vc  = 

C= 

Ua,Va  = 
a  = 
s  = 

Dx,Dy,Dz  = 

E  = 


depth 

x,y  velocity  components  at  the  water  surface 

water-surface  elevation 

x,y  velocity  at  the  bed 

bed  elevation 

salinity 

diffusion  coefficient  for  salt 
eddy  viscosity  components 


The  continuity  equation 


du  8v  dw 

—  +  — +  —  =  0 

dx  dy  dz 


(A6) 


is  solved  as  a  second  part  of  each  solution  step.  Equation  A6  is  converted  to  an 
appropriate  boundary  value  problem  through  differentiation  with  respect  to  z. 
After  rearrangement  it  takes  the  form 


d2  w  _  d  f  du  dv' 
dz2  dz\dx  dy , 


(A7) 


subject  to  boundary  conditions  specified  for  the  water  surface  and  the  bed. 


dC  d£  8h 

wr  =  ur  — 1 -  +  vr  — —  +  —  at  the  water  suface 

^  ^  dx  4  dy  dt 


(A8) 


and 


da  da  ,  ,  . 

wa  -  ua  ~  +  va  ~z  at  tae  bed 
dx  dy 


(A9) 


Note  that  in  these  equations  the  values  of  u  and  v  will  be  known  at  all  locations 
from  the  previous  part  of  the  solution  step.  Values  of  w  in  this  solution  are  used 
in  the  next  iteration  for  u,v,h,  and  s. 

The  geometric  system  varies  with  time;  i.e.,  the  water  depth  h  varies  during 
the  simulation.  In  order  to  develop  an  Eulerian  form  for  the  solution,  it  is  de¬ 
sirable  to  transform  this  system  to  one  that  can  be  described  with  a  constant 
geometric  structure.  Early  development  of  the  model  (King  1 982)  used  a  ex- 
transformation  in  which  the  bed  and  the  water  surface  are  transformed  to  con¬ 
stants.  In  a  later  analysis  of  this  method,  King  (1985)  pointed  out  that  at  locations 
where  a  sharp  break  in  bottom  profile  occurs,  the  transformation  is  not  unique 
and  momentum  in  the  component  directions  may  not  be  correctly  preserved.  An 
alternative  transformation  that  preserves  the  bottom  profile  as  defined,  but 


Appendix  A  Hydrodynamic  Code 


transforms  the  water  surface  to  a  constant  elevation  is  now  used  ( zv 
transformation). 

This  transformation  is  defined  by: 


X*=X 

(A10) 

II 

> 

(All) 

v  ,/  Jb-a) 

zv  =  a  +  (z-a) 

h 

(A12) 

where  b  is  the  fixed  vertical  location  to  which  the  water  surface  will  be  trans¬ 
formed.  Equations  A1-A6  and  A7-A9  then  incorporate  the  transformation  (A10- 
A12). 

Another  advantage  of  this  transformation  is  that  it  produces  zv  =  constant 
lines  that  are  close  to  horizontal,  i.e.,  z  =  constant  lines.  This  results  in  less 
fictitious  density-driven  currents  near  bed  profile  breaks  (Stelling  and  van  Kester 
1993).  Since  stratification-related  phenomena  are  usually  nearly  horizontal,  it  is 
important  that  the  transformation  leave  constant  surfaces  that  are  nearly 
horizontal.  Considering  the  pressure  gradient  (due  to  the  density  gradient)  in  this 
transformation  produces 


dP__dP_  dP  8ZV 
dx  dx  dz  dx 


(A13) 


In  a  strongly  stratified  stagnant  system  this  pressure  gradient  should  be  zero. 
However,  note  that  Equation  A1 3  in  the  transformed  system  is  dependent  upon 
two  terms  (each  of  which  could  be  large)  to  cancel  each  other.  This  could  cause 
artificial  currents  due  to  truncation  and  round-off  error.  A  transformation  in 
which  Qz v  I  dx  *  0  ,  i.e.,  zv  ~  z ,  will  reduce  this  problem.  Figure  A2  shows  an 
example  for  a  case  in  which  a  40-ft-deep  channel  passes  through  an  8-ft-deep  bay. 
Here  b  is  chosen  to  be  an  elevation  of  0  and  Cis  2  ft.  Near  the  break  in  the  bed 
profile  8zv/8x  is  fairly  small,  or  zv  surfaces  are  nearly  horizontal.  Contrast 
this  with  the  a  transformation  in  Figure  A3.  The  a  =  constant  surfaces  are  far 
from  horizontal  along  the  channel  side  slopes.  The  truncation  and  round-off 
errors  tend  to  drive  fictitious  currents  that  cause  the  denser  salt  water  to  leave  the 
channel.  The  transformation  results  in 


- f=o(C-b ) 


(A14) 


A4 
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Gteviikw,  Z 


whereas  the  a  transformation  is 


ir=0(h> 

ox 


which  is  much  larger. 


Figure  A2.  Lines  of  constant  z  near  a  significant  grade  change 
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Figure  A3.  Lines  of  constant  a  near  a  significant  grade  change 


The  Galerkin  finite  element  approximation  of  Equations  A1-A4  and  A7  uses 
a  quadratic  approximation  for  u,v,w,  and  s  and  linear  for  h  and  P.  The  non¬ 
linearity  is  addressed  by  Newton-Raphson  iteration  at  each  time-step.  Generally 
the  iteration  process  is  split  into  calculation  of  Equations  A1-A3,  then  A7, 
followed  by  A4.  This  sequence  is  repeated  until  sufficient  convergence  is 
reached. 
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